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SUMMARY 


The  satellite  NOAA-B  (1980-43A)  was  launched  in  May  1980  into  an 
orbit  with  perigee  height  near  260  km  and  apogee  height  1440  km,  at  an 
inclination  of  92.2°.  The  lifetime  was  11  months.  The  orbit  has  been 
determined  at  40  epochs  between  October  1980  and  May  1981  from  about  3000 
radar  and  optical  observations.  The  average  orbital  accuracy,  radial  and 
cross-track,  was  about  100  m,  with  rather  better  accuracy  in  the  final 
14  days. 

The  variation  of  orbital  inclination  has  been  analysed  to  determine 
two  good  values  of  atmospheric  notation  rate,  namely  1.10  *  0.10  rev/day 
at  300  km  (average  local  time)  and  1.15  t  0.06  rev/day  at  225  km  (evening). 

The  effect  of  atmospheric  rotation  on  the  precession  of  the  orbital 
plane  of  an  actual  satellite  has  never  previously  been  detected;  it  is 
clearly  apparent  for  1980-43A  in  its  last  days,  and  conforms  to  the 
expected  theoretical  change. 

The  variation  of  perigee  height  has  been  analysed  to  determine  ten 
values  of  atmospheric  density  scale  height,  for  heights  of  280-370  km. 

These  values,  accurate  to  about  3Z,  exceed  by  15Z  the  values  indicated  by 
the  COSPAR  International  Reference  Atmosphere .  Solar  activity  was  higher 
in  the  years  1980-81  than  at  any  time  since  early  1958,  and  it  appears 
that  the  CIRA  model  underestimates  the  density  and  density  scale  height  at 
high  levels  of  solar  activity. 
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1  INTRODUCTION 

The  meteorological  satellite  NOAA-B,  designated  1980-43A,  was  launched  on 
29  May  1980  but  failed  to  achieve  its  intended  orbit  and  decayed  within  a  year, 
on  3  May  198 1 1 .  The  satellite  was  a  rectangular  box,  3.71  m  long,  1.88  m  in 
diameter  and  just  over  1400  kg  in  weight,  of  which  the  fuel  supply  comprised 
almost  half.  Its  initial  perigee  and  apogee  heights  were  264  km  and  1445  km, 
and  the  orbital  inclination  was  92.23®. 

The  orbit  has  been  determined  at  40  epochs  from  about  3000  radar  and 

optical  observations  using  the  RAE  orbit  refinement  program  PROP  in  the  PROP  6 
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version  ;  orbits  were  determined  for  each  of  the  14  days  prior  to  decay  when  a 
large  number  of  observations  by  the  assigned  and  contributing  sensors  of  the 
North  American  Air  Defense  Command  (NORAD)  were  made  available  to  us. 

>The  main  purpose  of  this  Report  was  to  evaluate  the  atmospheric  rotation 
rate  from  changes  in  the  satellite's  orbital  Inclination.  This  was  successfully 
accomplished  and  it  was  also  possible  to  make  two  other  useful  analyses:  the 
right  ascension  of  the  node  was  analysed  to  demonstrate  for  the  first  time  the 
effects  of  atmospheric  rotation  on  this  orbital  element;  also,  since  the  solar 
activity  was  at  a  higher  level  than  at  any  time  since  1958,  the  values  of  scale 
height  derived  from  the  decrease  in  perigee  height  proved  to  be  of  unexpected 
significance.. 

\ 

2  ORBIT  DETERMINATION 

2. 1  Observations 

The  orbit  of  1980-43A  has  been  determined  at  40  epochs  from  12  October 
1980  until  its  decay.  The  observations  accepted  in  the  final  orbits  numbered 
2789,  consisting  of  over  1250  US  Navy  Navspasur  observations  supplied  by  the 

Naval  Research  Laboratory,  nearly  1100  NORAD  observations,  over  250  British 
radar,  168  visual  observations  from  volunteer  observers  reporting  to  the 
University  of  Aston  and  12  from  the  kinetheodollte  at  the  South  African 
Astronomical  Observatory. 

2.2  Observational  accuracy 

The  observational  accuracy  of  the  US  Navy  and  optical  observing  stations 
with  4  or  more  observations  accepted  is  indicated  in  Table  1  in  the  form  of  ras 
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residuals,  obtained  using  the  RAE  computer  program  ORES  . 
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Table  1 


Residuals  for  selected  stations 


Station 

Number 

of 

observations 

Rms  residuals 

Range 

km 

Minutes  of  arc 

RA 

Dec 

Total 

1 

US  Navy 

187 

2.3 

2.6 

3.4 

2 

US  Navy 

165 

3.0 

2.9 

4.2 

3 

US  Navy 

185 

2.8 

3.0 

4.1 

4 

US  Navy 

172 

2.8 

2.7 

3.8 

5 

US  Navy 

193 

2.3 

2.7 

3.5 

6 

US  Navy 

196 

3.2 

2.9 

4.3 

29 

US  Navy 

160 

0.8 

0.4* 

0.9* 

2125 

Street 

4 

1.7 

1.9 

2.5 

2265 

Farnhara 

5 

3.9 

4.0 

5.6 

2414 

Bournemouth 

60 

3.0 

4.2 

5.2 

2418 

Sunni ngdale 

4 

3.3 

2.5 

4.1 

2420 

Willowbrae 

78 

2.3 

2.0 

3.0 

2431 

Copthorne 

4 

0.7 

0.7 

1.0 

2572 

Wittenheim 

4 

3.2 

5.0 

6.0 

2577 

Cape  Kinetheodolite 

12 

0.58 

0.65 

0.87 

♦Geocentric 


The  residuals  represent  the  sum  of  the  observational  errors  and  orbital 
errors.  Since  the  drag  was  severe  and  variable,  the  orbital  errors  may  occa¬ 
sionally  be  Important.  Certainly  the  residuals  for  both  the  US  Navy  and  optical 
observations  are  slightly  larger  than  usual,  and  this  seems  the  most  likely 
explanation. 

The  satellite's  orbit  was  nearly  sun-synchronous  and  this  led  to  long 
intervals  of  invisibility,  thus  reducing  the  number  of  visual  observations.  As 
in  many  previous  orbit  determinations,  the  observations  from  Willowbrae 
(R.D.  Eberst)  and  Bournemouth  (D.J.  Hopkins)  made  up  a  very  large  proportion  of 
the  visual  observations. 

2.3  Orbital  elements  and  accuracy 

The  computed  sets  of  orbital  elements  with  their  standard  deviations  at 
each  of  the  40  epochs  are  listed  in  Table  2.  The  epoch  for  each  orbit  is  at 
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00  hours  on  the  day  indicated  and  the  PROP  program  fits  the  mean  anomaly  M  by  a 
polynomial  of  the  form 

M  =  MQ  +  Mjt  +  M2t2  +  Hjt3  +  M4t**  +  Mjt5  ,  (1) 

where  t  is  the  time  measured  from  epoch  and  the  number  of  M-coef f Icients  used 
depends  on  the  drag  experienced  by  the  satellite  at  around  that  epoch.  For 
1980-43A,  the  coefficients  Mq  -  were  sufficient  for  just  3  of  the  40 
orbits,  15  orbits  required  M2  to  be  added,  16  orbits  needed  M4  and  6  orbits 
required  the  full  complement  of  Mq-  .  The  measure  of  fit,  shown  by  £  , 
varied  between  0.30  and  1.36  with  an  average  value  of  0.67,  indicating  that  the 
observations  were  fitted  well  on  most  of  the  orbits. 

The  values  of  inclination,  i  ,  plotted  in  Fig  1,  had  standard  deviations 
between  0.0004°  and  0.0024°,  the  average  being  0.0012°,  corresponding  to  about 
140  m  in  distance.  For  the  final  14  dally  orbits,  the  sd  varied  from  0.0004°  to 
0.0016°  (the  latter  on  the  date  of  decay),  with  an  average  of  0.0007°.  The 
right  ascension  of  the  node,  ft  ,  had  standard  deviations  202  lower  than  i  ,  on 
average. 

For  eccentricity,  e  ,  the  sd  ranged  between  4  x  10  6  and  3  x  10  5,  with 

an  rms  of  1.4  x  10  5  or  100  m  in  distance,  while  the  rms  value  for  the  final  14 

orbits  was  1  x  10  5  (70  m).  The  values  of  a(l  -  e)  from  Table  2  are  plotted 

in  Fig  2.  They  show  the  expected  oscillation  due  to  odd  zonal  harmonics,  of 
amplitude  about  10  km,  and  the  general  decrease  due  to  drag  as  the  lifetime 
progresses. 

Fig  3  shows  the  values  of  M2  ,  which  is  proportional  to  the  drag  acting 
on  the  satellite.  There  are  many  irregular  variations,  which  reflect  the 
response  of  air  density  to  the  high  and  variable  solar  activity  prevailing  in 
the  solar-maximum  period  of  1980-81.  Underlying  the  irregular  changes  is  a 
smoother  variation  influenced  by  the  oscillations  in  perigee  height,  particu¬ 
larly  the  final  steep  decrease  after  MJD  44680.  Changes  in  M2  may  also  have 
occurred  because  of  changes  in  the  satellite's  cross-sectional  area. 

3  ZONAL  WINDS  FROM  ANALYSIS  OF  INCLINATION 
3.1  Introduction 

The  decrease  of  inclination  during  the  satellite's  life,  visible  in  Fig  1, 
is  caused  mainly  by  the  aerodynamic  forces  normal  to  the  orbital  plane,  created 
by  atmospheric  rotation.  If  other  perturbations  are  removed,  the  decrease  in  1 
can  be  analysed  to  determine  the  atmospheric  rotation  rate,  which  gives  the 
zonal  (west-to-east)  wind. 
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3.2  Removal  of  perturbations 

The  lunisolar  and  zonal  harmonic  perturbations  have  been  calculated  using 
the  computer  program  PROD*4,  and  the  perturbations  due  to  the  ^2,2  tesseral 
harmonic  are  given  with  the  PROP  printout.  Removal  of  these  sources  of 
perturbation  gives  the  values  of  Inclination  plotted  in  Fig  4. 

There  are  further  perturbations  in  inclination  due  to  earth  and  ocean 
tides,  and  the  change  in  reference  system  produced  by  the  precession  of  the 
Earth's  axis.  The  magnitude  of  these  perturbations  for  polar  orbits  is  indi¬ 
cated  in  Ref  5:  the  largest  earth-tide  effect  is  that  of  lunar  origin,  which  can 
change  the  inclination  by  as  much  as  0.0002“  per  100  days;  this  effect  is 
negligible  here,  since  the  smallest  standard  deviation  in  inclination,  in  the 
first  190  days  is  0.0007°.  The  correction  to  i  required  as  a  result  of  pre- 

2 

cession  is  given  approximately  by 

—4 

AiDOK..  -  8,7  —  (cos  a  -  cos  ».)  degrees  (2) 

rWiL  ^  1 

where  flj  and  SI  are  the  final  and  current  values  of  Si  ,  and  SI  is  in 
degrees/day.  Here  ft  =  0.28  deg/day,  so  the  numerical  factor  in  (2)  is 

3.1  x  10  3.  During  the  last  14  days  of  the  satellite’s  life 
and  is  negligible.  During  the  first  190  days  the  maximum  value  of 
|  cos  0  -  cos  S^  |  is  0.17,  so  that  AlpREC  <  0.0005°,  and  this  is  also  negli¬ 
gible,  since  the  average  sd  of  the  values  of  i  in  the  first  190  days  is 
0.0012°. 

Perturbations  due  to  meridional  winds  are  negligible  because  the  inclina¬ 
tion  is  so  close  to  90°. 

Perturbations  due  to  resonance  with  the  Earth's  gravitational  field  may  be 
significant,  and  will  be  discussed  in  section  3.3.2. 

3.3  Analysis 

3.3.1  Between  12  October  1980  and  19  February  1981 

The  theoretical  variation  of  1  with  time  for  a  series  of  values  of 
atmospheric  rotation  rate  A  ,  was  obtained  using  the  ROTATM  computer  program, 
and  the  value  of  A  ■  1.10  rev/day  best  fits  the  observed  points  on  orbits  1-19, 
between  12  October  1980  and  19  February  1981,  as  shown  in  Fig  4. 

As  the  scale  at  the  top  indicates,  the  local  time  at  perigee  is  initially  g 

at  08-09  h,  then  jumps  by  10  hours  as  perigee  passes  near  the  south  pole  on  * 

3 

3  November,  and  remains  near  17-18  h  until  25  December,  when  perigee  passes  near  ’ 

i 

* 
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the  north  pole  and  the  local  time  shifts  to  03-04  h  until  mid-February  1981.  It 
is  usually  helpful  to  classify  the  values  of  A  obtained  from  orbit  analysis 
according  to  their  local  time  in  one  of  three  categories,  namely,  morning 
(06-12  h),  evening  (18-24  h)  or  average.  This  first  value  of  A  obtained 
here  is  obviously  not  biased  towards  either  morning  or  evening,  and  falls  into 
the  'average'  category.  Since  the  perigee  moves  from  high  southerly  latitudes 
across  the  equator  to  the  north  pole  and  back  again,  between  October  and 
February,  there  is  no  significant  seasonal  bias. 

As  all  the  19  points  are  less  than  0.0024®  from  the  curve,  it  is  reasonable 
to  assign  a  standard  deviation  of  0.0012®  to  the  curve,  and  so  the  total 
decrease  in  inclination,  namely  0.019®,  is  likely  to  be  measured  with  an  error 
of  about  /2  x  0.0012®,  that  is  0.0017®  or  9%.  Thus  the  appropriate  accuracy  for 
A  is  1.10  ±  0.10. 


The  average  perigee  height  during  this  time  is  260  km,  and  the  height  at 
which  the  value  of  A  applies  may  be  taken  as  iH  above  this,  where  the 
density  scale  height  H  is  about  55  km.  So,  the  final  result  is: 

A  “  1.10  ±  0.10  for  300  km  height,  average  local  time  and  average  season. 


3.3.2  Between  19  February  and  2  April  1981 


The  satellite  experienced  l5th-order  resonance  with  the  Earth's  gravita¬ 
tional  field  on  4  March  1981.  Because  of  the  high  drag  at  this  time,  the 
passage  through  resonance  was  rapid,  but  the  resulting  perturbation  still  seemed 
likely  to  be  appreciable. 


So  an  extra  orbit  was  determined  for  5  March  1981  and  further  orbits  were 


determined  by  PROP  for  epochs  two  days  before  and  after  this  orbit  and  before 
and  after  orbits  20,  21  and  22,  and  three  days  after  orbit  19,  to  give  a  total 
of  14  orbits  during  the  period  when  the  inclination  was  likely  to  be  affected  by 
resonance.  The  change  in  Inclination  due  to  resonance  was  then  calculated  using 
the  THROE  computer  program  with  values  for  lumped  harmonics  taken  from  Ref  7, 
namely: 


- 

io9c;^  - 

-16.8 

io9S°^  - 

-5.4 

.  ‘  <3> 

* 

- 

-52 

1°  ST5  “ 

-42 

j 

The  effects  of  the  (C.S)1’0 

coefficients 

are  negligible 

at  this  inclination. 
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The  resulting  variation  of  i  due  to  resonance  is  shown  in  Fig  5,  where 
the  initial  value  of  1  is  arbitrary  and  only  the  change  in  i  is  of  interest. 
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The  pre-resonance  oscillation  in  Fig  5  is  centred  on  a  value  close  to  92.1931° 
and  the  post-resonance  oscillation  is  centred  on  a  value  close  to  92.1961°,  so 
that  the  resonance  produces  an  increase  of  0.0030°  in  i  • 

The  overall  variation  in  i  ,  shown  as  a  broken  line  in  Fig  4,  has  been 
calculated  by  adding  (i  -  92.1931)  from  Fig  5  to  the  values  given  by  ROTATM,  the 
value  of  A  in  ROTATM  (0.825)  being  chosen  to  bring  the  curve  to  the  correct 
value  for  the  start  of  the  final  section,  from  2  April  1981  onwards,  namely 
i  -  92.1787. 

The  value  of  0.825  required  for  A  is  not  accurate  enough  to  qualify  as  a 
'determination'  of  A  ,  since  it  has  a  standard  deviation  of  about  0.15.  As 
the  local  time  at  perigee  is  12-13  h,  there  is  a  bias  towards  morning  (06-12  h) 
when,  according  to  Ref  6,  a  value  of  A  near  0.85  would  be  expected,  for  height 
280  km  and  average  season. 

In  view  of  the  high  drag,  the  change  in  inclination  due  to  the  15th-order 
harmonics  in  the  geopotential  is  surprisingly  large,  and  it  is  clear  from  the 
form  of  Fig  5  that  nearly  the  largest  possible  increase  has  happened  to  occur. 
The  perturbation  at  the  31:2  resonance  (13  April)  should,  under  the  usual 
assumption  of  proportionality  to  ^  ,  be  less  than  one  tenth  of  that  at 
l5th-order  resonance,  that  is,  less  than  0.0003“  and  hence  negligible.  The  29:2 
resonance  occurs  in  October  1980  between  orbit  2  and  orbit  3,  and  on  the  same 
assumptions  could  produce  a  maximum  change  in  inclination  of  0.0007°.  This  can¬ 
not  be  regarded  as  important,  because  the  accuracy  of  the  curve  at  this  time  has 
been  estimated  as  0.0012°;  but  if  a  net  decrease  of  0.0007°  did  occur  at  29:2 
resonance,  the  first  two  points  in  Fig  4  would  fit  better.  Good  values  for 
lumped  coefficients  of  order  29  are  not  available,  so  the  change  cannot  be 
reliably  calculated. 

3.3.3  Between  2  April  1981  and  decay  (3  May  1961) 

Fig  6  shows  the  values  of  Inclination  and  the  fitted  A  curve  in  the 
final  phase  of  the  satellite's  life.  (The  curve  is  also  given  in  Fig  4  on  a 
smaller  scale.)  The  value  of  A  required  is  1.15.  The  fitting  is  reasonably 
good,  and  there  is  no  obvious  break-point  such  that  a  better  fitting  could  be 
achieved  by  using  two  different  values  of  A  .  The  standard  deviation  in  A  is 
difficult  to  estimate,  because  the  values  of  Inclination  at  the  beginning  and 
end  are  less  accurate  than  the  Intermediate  values.  If  we  assign  an  error  of 
0.002°  at  the  beginning  and  0.001°  at  2  May,  the  sd  in  A  is  0.06  since  the 
total  decrease  in  inclination  is  0.042°.  Alternatively,  an  error  of  0.001° 
might  be  assigned  between  20  April  and  2  May,  where  the  decrease  in  inclination 
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is  0.0275s;  this  also  gives  the  sd  of  A  as  0.06.  So  we  quote 
A  -  1.15  ±  0.06. 


Fig  4  shows  that  the  local  time  at  perigee  is  between  22  and  24  h,  so  the 
value  of  A  is  in  the  'evening'  category  of  Ref  6.  Since  perigee  moves  from 
latitude  88s  north  on  2  April  to  45s  south  on  3  May,  the  season  is  'average'. 
The  perigee  height  decreases  from  220  km  on  20  April  to  182  km  on  2  May  (the 
time  interval  where  there  are  accurate  daily  values),  so  the  average  perigee 
height  is  about  200  km,  and  the  height  at  which  A  applies  is  about  225  km. 
Thus  we  have: 


A  =  1.15  ±  0.06  for  225  km  height,  evening  and  average  season. 
4  ANALYSIS  OF  RIGHT  ASCENSION 


4.1  Introduction 


The  aerodynamic  force  produced  by  atmospheric  rotation  causes  pertur¬ 
bations  in  the  right  ascension  of  the  ascending  node  ft  ,  as  well  as  in  the 
inclination.  However,  these  atmospheric  perturbations  in  0  are  usually  very 
small  and  they  cancel  out  during  each  half  cycle  of  the  argument  of  perigee  m  . 
Consequently,  the  atmospheric  effect  on  right  ascension  has  never  been  conclus¬ 
ively  demonstrated.  The  orbit  of  1980-43A  in  its  last  few  days  offers  some  hope 
of  successful  analysis  of  12  ,  and  this  possibility  is  pursued  below. 


4.2 


The  effect  of  atmospheric  rotation  on  11 

8 

The  change  in  12  due  to  an  atmosphere  rotating  at  a  rate  A  is  given  by 


AO 


A  sin  2(d 
6/F 


©I* 


cos  2(0  +  0(c2,e2) 


AT  ,  (4) 


where  I  is  the  Bessel  function  of  the  first  kind  and  imaginary  argument,  of 
n  _ 

argument  z  “  ae/H  ;  /F  is  a  factor  which  has  a  value  1.00  for  1980-43A;  and 
AT  is  the  change  in  orbital  period  expressed  in  days.  The  parameter  c 
expresses  the  effects  of  atmospheric  oblateness  e'  ,  being  given  by 
c  “  {e'a(l  -  e)  sin2  i}/2H. 


For  most  satellites  n  changes  rapidly,  at  a  rate  of  several  degrees  per 
day,  as  a  result  of  the  gravitational  forces  produced  by  the  departure  of  the 
geopotential  from  spherical  symmetry.  For  1980-43A,  however,  which  is  a  near- 
polar  orbit,  the  rate  is  only  0.3s  per  day  and  there  is  hope  of  removing  this 
perturbation  without  significant  error. 

If  a  measurable  change  in  12  is  to  build  up,  it  is  obvious  from  equation 
(4)  that  thi  a  must  '  -  large  change  AT  in  the  orbital  period,  while  sin  2w 


-  ■  r 
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retains  the  same  sign.  For  1980-43A,  the  most  promising  time  interval  is  the 
last  8  days  of  the  life,  where  the  largest  change  AT  occurs  and  sin  2us  is 
consistently  negative  (see  Table  2). 

4.3  Procedure  for  removal  of  perturbations 

The  terrestrial  gravitational  perturbations  in  0  were  calculated  at 
daily  intervals,  with  a  re-start  each  day,  using  the  computer  program  PROD  and 
harmonics  up  to  degree  20.  Since  PROD  does  not  include  the  effects  of  drag, 
this  procedure  gives  8  daily  values  of  0  due  to  terrestrial  gravity  at  dates 
April  26.0,  27.0  ...  May  3.0.  The  first  7  of  these  were  integrated  with  the 
Gregory  central-difference  formula  for  integration  interval  w  , 

w 

I  hdt  «  wfy  "TT  +T5o  +  •••)  ,  (5> 


where  ft 


* 


+  V’ 


AUn 


a  - 

w  0 1 


etc. 


620. 


KA2S2 


+  A2ii0)  ,  etc. 


Here  w  =  1  and  the  terra  is  negligible.  The  running  total  of  these 

integrals,  which  we  denote  by  A flg  ,  gives  the  total  perturbation  to  be 
subtracted  from  the  PROP  value  of  «  (Table  2)  in  order  to  clear  terrestrial 
gravitational  perturbations. 

Equation  (5)  cannot  be  applied  between  May  2.0  and  3.0,  because  S2  -*■  * 
early  on  May  3  and  the  central  difference  value  becomes  unreliable.  So  an 

analytical  integration  is  necesary.  If  suffix  0  denotes  values  at  May  2.0, 

0 

the  orbital  period  T  at  a  later  time  is  given  by  theory  as 


T  .  3H  .  z0Il(z0) 

Tq  “  1  “  2a  in  zIL(z) 


(6) 


Between  May  2.0  and  3.0,  the  value  of  z  decreases  from  about  1.4  to  about  1.0, 

a  region  in  which  I.(z)  is  approximately  proportional  to  z^"*  .  Also  the 
8  2  2 

theory  indicates  that  z  /zq  =  1  -  At  ,  where  t  is  the  time  after  May  2.0, 
and  1/A  is  approximately  the  lifetime  after  May  2.0.  Thus  equation  (6)  gives 

i;  ■  1  -  If  (r)  ■  in  U  -  At)1-05  .  (7) 

”3.5 

The  rate  of  change  of  0  due  to  terrestrial  gravity  is  proportional  to  a 
-7/3 

and  hence  to  T  •  Therefore 


CJ 

o 

o 


Wf 
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K(T/T0)_7/3  =  K{  1  -  3^-  In  (1  -  At)}  , 


1  I 


A  - 


(8) 


where  K  Is  constant.  The  PROD  output  shows  that  h  *  0.3311  deg/day  at 
May  2.0  and  0.3393  deg/day  at  May  3.0.  Substituting  these  values  in  equation 
(8),  and  taking  llH/3.a  as  0.016,  we  find  A  =  0.788.  Hence  the  perturbation  in 
0  due  to  terrestrial  gravity  from  May  2.0  to  3.0  is 

1 

I  0.3311 { 1  -  0.016  In  (1  -  0.788t)}dt  =  0.3342°  .  (9) 

0 


The  lunlsolar  perturbation  to  0  was  calculated  by  re-running  the  PROD 
program  (with  daily  restart)  with  lunisolar  perturbations  included,  and  then 
subtracting  the  previous  values  with  lunisolar  terms  absent.  These  accumulated 
lunisolar  perturbations,  say,  were  also  subtracted  from  the  original  PROP 

values,  SI 


PROP' 


Thus  we  arrive  at  a  value  S20  free  of  terrestrial  and 


lunisolar  gravitational  perturbations  and  given  by 


nPR0P  ”  AiJG  MJLS 


4.4  Results 

The  values  of  are  plotted  in  Fig  7  together  with  a  curve  giving  the 

theoretical  change  due  to  an  atraospnere  rotating  at  1.15  rev/day,  the  value 
appropriate  for  the  last  7  days  of  the  satellite's  life  (Fig  6).  Although  the 
total  expected  change  in  SI  due  to  atmospheric  rotation  during  the  7  days  is 
only  0.0028°,  the  values  of  S^  and  the  curve  are  in  good  agreement,  none  of 
the  points  being  more  than  2  sd  from  the  curve. 

Thus,  the  change  in  S2  due  to  atmospheric  rotation  is  consistent  with  an 
atmospheric  rotation  rate  of  A  •  1.15  .  The  values  are  not  accurate  enough  to 
allow  an  independent  determination  of  A  of  useful  accuracy,  but  the  agreement 
between  theory  and  observation  is  close  enough  to  justify  the  conclusion  that 
the  observational  standard  deviations  are  realistic  and  that  the  theory  -  which 
has  not  been  previously  tested  against  observation  -  is  basically  correct. 

5  ANALYSIS  OF  PERIGEE  HEIGHT 

5.1  Determining  density  scale  height:  procedures  and  results 

Fig  2  shows  that  the  perigee  distance  a(l  -  e)  suffers  considerable 
purturbations  due  to  the  odd  zonal  harmonics  in  the  geopotential,  and  smaller 
lunisolar  gravitational  perturbations.  These  have  been  removed  using  the  PROD 
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computer  program  with  daily  integration  steps  and  restarts  at  intervals  of  about 
50  days,  to  give  a  perigee  distance  freed  of  gravitational  perturbations.  On 
subtracting  the  mean  Earth  radius,  6367  km,  we  obtain  a  perigee  'height*  Q  free 
of  gravitational  perturbations,  plotted  in  Fig  8. 

The  values  of  Q  in  Fig  8  show  a  continual  decrease  due  to  the  effects  of 

air  drag,  and,  if  z  >  3,  the  rate  of  decrease  of  Q  in  an  atmosphere  of 

g 

ellipticity  e'  may  be  expressed  as 


Hl“l  /.  _  ,  .  H  2e  ' 

SM^e  y  6  4ae  e 


sin  i 


2u) 


(11) 


If  H  denotes  the  value  of  the  density  scale  height  at  perigee,  the  quantity 

Hj  in  equation  (11)  is  the  value  of  density  scale  height  at  a  height  above 

perigee. 

Values  of  Q  for  use  in  equation  (11)  are  found  in  the  form  AQ/At  , 
where  the  time  interval  At  is  chosen  to  be  long  enough  to  ensure  that  the 
change  AQ  in  Q  is  large  enough  to  be  determined  with  a  formal  accuracy 
better  than  3%.  Since  the  standard  deviations  in  the  values  of  a(l  -  e)  are 
of  order  100  m,  the  value  of  AQ  should  be  at  least  3  km.  The  values  of  AQ 
are  obtained  from  the  mean  curve  in  Fig  8,  rather  than  the  individual  points. 

The  values  of  ,  Mj ,  e  and  cos  2 w  for  use  in  equation  (11)  are  of  course 
also  averaged  over  the  time  interval  At  .  This  averaging  should  not  lead  to 
errors  of  more  than  1%:  the  e'  term  is  normally  small  and  this  much  reduces 
the  effect  of  errors  in  averaging  cos  2w  . 

The  use  of  equation  (11)  for  calculating  H.  from  the  observational 
•  •  *■ 

values  of  Q,  ,  e,  etc  has  the  advantage  that  it  does  not  require  any 

'model  atmosphere'  to  be  specified.  The  only  assumption  is  that  the  atmosphere 
is  oblate,  with  the  ellipticity  e'  taken  as  0.00335. 

The  values  of  Hj  obtained  from  equation  (11)  are  shown  in  histogram  form 
in  Fig  9  as  unbroken  horizontal  lines.  For  the  last  two  values  of  in  Fig  9, 

equation  (11)  is  nearing  the  limits  of  its  validity  because  z  **  ae/Hp 
decreases  to  near  3,  and  a  correction  is  needed.  When  z  =*  3,  the  assumption  of 
an  oblate  atmosphere  symmetrical  about  the  Earth's  axis  may  be  over-simple,  and 
there  is  the  possibility  of  a  significant  contribution  to  equation  (11)  from  the 
asymmetry  produced  by  the  day-to-night  variation  in  air  density.  If  f  denotes 
the  ratio  of  maximum  daytime  density  to  minimum  night-time  density, 

F  “  (f  -  l)/(f  +  1),  and  $p  is  the  Sun-perigee  angle  (or,  strictly,  the 
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geocentric  angle  between  perigee  and  the  daytime  density  maximum),  the  effect 
of  the  day-to-night  density  variation  may  be  expressed  in  terms  of  the  para¬ 
meter  u  defined  by  u  =  (F  cos  *  )/(l  +  F  cos  6  )  .  To  take  account  of  this 

P  P 

effect,  it  is  necessary  to  add  an  extra  terra, 


j 

j 

o 


within  the  brackets  in  equation  (11).  This  term  was  found  to  be  appreciable  for 
the  last  two  values  of  in  Fig  9,  and  was  taken  into  account  in  their  calcu¬ 

lation.  It  is  not  significant  for  the  previous  value,  for  which  <J>  *  90°  ,  nor 
for  earlier  values,  when  z  is  much  larger.  With  this  correction  it  is  hoped 
that  an  accuracy  of  3%  in  has  been  maintained. 

The  perigee  height  y^  is  obtained  by  subtracting  the  local  Earth  radius 

from  a( 1  -  e),  and  the  height  y^  at  which  applies  (shown  in  Fig  9)  is 

taken  as  y  +  4h  ,  where  H  is  determined  as  described  in  the  next  section. 

P  P  P 

(The  small  correction  to  yp  ,  dependent  on  ui  ,  has  not  been  applied  because 
greater  errors  in  yj  arise  from  errors  in  Hp.) 

The  evaluation  of  H  has  not  been  continued  beyond  20  April,  because  the 

increased  importance  of  the  oblateness  and  day-to-night  terms  would  prevent  the 

3%  accuracy  being  maintained.  (To  be  specific,  the  e'  term  in  equation  (11) 
has  a  value  of  0.13  for  the  last  point  in  Fig  9,  and  would  have  a  value  of  0.3 
for  the  next  (uncalculated)  point;  and  t'  may  have  errors  of  10-20X  in  the 
geomagnetlcally  disturbed  week  of  21-27  April.)  Values  of  H  for  21  April- 
3  May  from  satellite  orbit  analysis  can  be  found  in  Ref  10. 

5.2  Comparisons  with  CIRA  values 

The  COSPAR  International  Reference  Atmosphere  1 972 1 1  has  been  widely  used 

as  a  standard  for  comparison  in  recent  years  and  has  proved  quite  reliable, 

apart  from  some  deficiencies  in  the  representation  of  the  semi-annual  and 
i  2 

geomagnetic  effects  .  Except  for  a  few  results  from  1958-9  at  only  one  height, 
CIRA  1972  was  based  on  the  mass  of  results  obtained  during  the  1960s,  when  the 
solar  activity  index  never  exceeded  180,  in  the  usual  units 

(10  22  w  m  2  Hz  1).  During  the  lifetime  of  1980-43A,  the  values  of  , 

which  are  shown  in  Fig  9,  averaged  about  220  with  three  week-long  excursions 
above  250.  Thus  1980-43A  provides  a  test  of  CIRA  1972  outside  its  normal  range. 

The  main  Tables  of  CIRA  give  values  of  density  and  density  scale  height 
for  heights  of  110  to  2000  km  for  a  series  of  values  of  exospheric  temperature 
Seven  parameters  are  required  in  order  to  calculate  T„: 
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(1)  the  3-monthly  averaged  value  of  7* 

(2)  the  current  value  of  Sl0>7,  given  In  Fig  9; 

(3)  the  current  value  of  the  geomagnetic  planetary  Index  Ap,  also  given 
in  Fig  9; 

(4)  the  semi-annual  variation,  specified  in  graphs  and  tables  in  C1RA  as 
a  function  of  day-of-the-year; 

(5)  -  (7)  the  local  time,  latitude  and  season  appropriate  to  perigee; 

these  have  been  calculated  from  the  orbital  elements. 

From  the  values  of  T—  thus  obtained  (they  range  between  1200  and  1420K), 

we  find  the  CIRA  value  of  density  scale  height  at  perigee  height  y^  . 

This  value  of  H  leads  to  a  value  of  y,  =  y  +  and  the  CIRA  value  of 

P  1  P  P  - 

at  height  y^  can  then  be  found.  The  CIRA  values  are  shown  as  a 

broken-line  histogram  at  the  bottom  of  Fig  9. 

It  is  obvious  from  Fig  9  that  the  values  of  obtained  from  analysis  of 

1980-43A  are  significantly  higher  than  the  CIRA  values;  it  is  also  clear  that 
the  discrepancy  becomes  smaller  as  time  goes  on,  presumably  as  a  result  of  the 
decrease  in  height  from  370  km  initially  to  282  km  at  the  last  point. 

For  the  first  five  points  (y^  £  339  km)  the  values  of  Hj  from  1980-43A 
are  on  average  18%  higher  than  the  CIRA  values,  a  margin  that  greatly  exceeds 
the  expected  3%  error  in  the  former.  For  the  other  five  points,  at  heights  of 
280-320  km,  the  difference  is  less,  12%  on  average,  but  still  highly  signifi¬ 
cant.  So  it  seems  that  the  values  of  density  scale  height  obtained  from 
CIRA  1972  for  heights  of  280-370  km  are  considerably  too  low  when  the  solar 
activity  index  Sjq>7  is  at  levels  near  220. 

It  could  be  argued  that,  since  the  CIRA  values  of  H  appear  to  be  too 
low,  the  value  of  H  used  in  calculating  y^  should  also  be  increased  by 
12-18%.  This  would  increase  y^  and  reduce  the  discrepancy,  which  would  become 
10-16%  rather  than  12-18%. 

The  significance  of  the  results  is  discussed  in  section  6.3 
6  DISCUSSION  OF  THE  RESULTS 
6. 1  Zonal  winds  from  analysis  of  inclination 

The  orbit  of  1980-43A  was  particularly  suitable  for  determining  zonal 
winds  from  analysis  of  inclination,  first  because  the  initial  period  was  high 
(100  minutes)  and  second  because  the  orbit  was  nearly  polar,  thus  leading  to  the 


ft 


15 


greatest  possible  change  in  inclination  (in  fact  nearly  0.09°)  and  also  ensuring 
that  any  uncertainty  created  by  the  effects  of  meridional  winds  was  negligible. 
There  was  only  one  resonance  (15th  order)  with  an  appreciable  effect,  and  that 
was  accurately  calculated  (Fig  5).  The  resulting  variation  of  inclination  with 
time  (Figs  4  and  6)  yields  two  values  of  atmospheric  rotation  rate  that  should 
be  completely  reliable,  namely: 

A  ”  1.10  ±  0.10  at  height  300  km,  average  local  time,  average  season; 

A  ■  1.15  ±  0.06  at  height  225  km,  evening  and  average  season. 

These  values,  which  imply  west-to-east  winds  of  40  ±  40  m/s  and  60  ±  24  m/s 
respectively,  have  already  been  used  in  a  recent  review  of  upper-atmosphere 
winds6. 

6.2  Analysis  of  right  ascension 

The  effect  of  atmospheric  rotation  on  right  ascension  has  long  been  known 
1 3 

on  theoretical  grounds  ,  and  Indeed  was  calculated  as  amounting  to  0.04°  for 
Sputnik  2.  This  was  very  much  smaller  than  the  change  in  right  ascension  due  to 
the  gravitational  effect  of  the  earth's  flattening,  which  was  more  than  500°,  so 
atmospheric  rotation  did  not  need  to  be  taken  into  account  when  using  Sputnik  2 
for  determining  the  even  zonal  harmonics  in  the  earth's  gravitational  field. 

The  effect  of  atmospheric  rotation  on  right  ascension  did  not  need  to  Vs  taken 
into  account  subsequently,  because  no  other  satellite  of  such  high  drag  as 
Sputnik  2  has  been  used  in  studies  of  gravitational  field. 

The  actual  change  in  right  ascension  produced  by  atmospheric  rotation  has 
not  previously  been  analysed,  because  equation  (4)  shows  that  a  successful 
analysis  could  only  be  made  if  several  diverse  conditions  were  satisfied: 

(a)  The  drag  must  be  very  severe,  so  that  a  large  change  in  period  AT 
could  build  up  quickly. 

(b)  The  value  of  sin  2u  has  to  remain  large  and  of  the  same  sign,  which 
usually  means  that  the  time  scale  should  not  be  more  than  a  week  or  two. 

(c)  The  orbit  must  not  be  too  nearly  circular  -  otherwise  Ij/Iq 
be  very  small. 

There  are  two  other  obvious  conditions: 

(d)  The  orbit  must  be  accurately  determined. 

(e)  The  perturbations  due  to  even  zonal  harmonics  mist  be  accurately 
removed  -  a  difficult  task,  since  they  can  amount  to  30°  or  more  unless 
the  orbit  is  near-polar. 
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All  of  these  five  criteria  are  unlikely  to  be  fulfilled  together,  but  they 
are  nearly  satisfied  by  1980-43A  during  its  last  week  in  orbit;  its  near-polar 
orbit  reduces  the  gravitational  effect  to  only  a  few  degrees  and  the  main 
weakness  is  the  low  eccentricity  -  a  higher  value  would  have  given  a  larger 
change  and  better  results* 

The  total  theoretical  change  in  right  ascension  for  1980-43A  during  its 
last  8  days  is  0.0028s  and,  as  Fig  7  shows,  theory  and  observation  are  in 
excellent  agreement,  thus  clearly  demonstrating  the  effect  for  the  first  time, 
confirming  the  theory  and  validating  the  standard  deviations  of  the  observa¬ 
tional  values. 

6.3  Density  scale  height  from  analysis  of  perigee  height 

It  is  not  possible  to  calculate  the  air  density  from  analysis  of 
1980-43A,  because  its  mass  and  cross-sectional  area  are  not  accurately  known  and 
the  cross  section  may  vary.  But  it  is  feasible  to  calculate  the  density  scale 
height  H  from  the  decrease  in  perigee  height  because  this  does  not  depend  on 
either  the  value  or  the  constancy  of  the  mass/area  ratio.  Previous  deter¬ 
minations  of  density  scale  height  from  satellite  orbits  (eg  Refs  9  and  14)  have 
given  values  which  are  on  average  within  42  of  those  obtained  from  the  COSPAR 
International  Reference  Atmosphere,  CIRA  1972. 

The  values  of  H  from  analysis  of  1980-43A,  which  should  be  accurate  to 

about  32,  are  much  higher  than  those  given  by  CIRA  1972  -  182  higher  on  average 

for  the  first  five  points  at  heights  above  335  km,  and  122  higher  on  average  for 

the  last  five  points  at  heights  of  280-320  km.  We  conclude  that  CIRA  1972 

seriously  underestimates  the  density  scale  height  at  these  high  levels  of  solar 
activity  (Sjq  i  *  220)  at  heights  of  280-370  km. 

This  result  is  not  too  surprising  because  CIRA  1972  is  known  to  under¬ 
estimate  the  effect  of  strong  transient  solar  disturbances.  Consistently  high 
levels  of  solar  activity  (Sjq  ^  >  200)  had  not  occurred  since  1958  and  CIRA 
1972  contains  very  little  data  from  this  era.  Since  the  density  scale  height 
is  proportional  to  the  temperature  divided  by  the  molecular  weight,  it  is  pro¬ 
bable  that  the  exospheric  temperature  T  is  higher  than  Indicated  by  CIRA 
1972.  (A  lower  molecular  weight  is  less  likely,  because  increased  temperature 
generally  goes  with  a  higher  molecular  weight). 

The  density  scale  height  H  is  inversely  proportional  to  the  height 
gradient  of  density  p  ,  being  defined  by  the  equation  dp/dy  ■  -p/H  ,  where  y 
is  height;  so  lower  values  of  H  imply  a  faster  decrease  in  p  as  height 
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increases,  therefore  it  is  to  be  expected  that  the  low  values  of  H  in 

CIRA  1972  will  lead  to  low  values  of  density  at  these  levels  of  solar  activity. 

This  expectation  is  confirmed,  though  not  conclusively,  by  the  value  of  density 

obtained  from  the  orbit  of  Sputnik  1  in  October  1957  when  the  monthly-averaged 

S  was  281.  The  value  of  density  from  Sputnik  1,  as  recently  revised15,  is 

10 « 7 

(1.7  ±  0.2)  x  10_l°  kg/m  at  247  km,  whereas  CIRA  1972  (with  T„  -  1500  K)  gives 
1.4  x  10  10  kg/m^.  An  increase  of  T^  to  1800  K  would  give  1.7  x  10  10. 

So  we  may  tentatively  conclude  that: 

(a)  the  density,  as  well  as  the  density  scale  height,  is  probably  too  low 
in  CIRA  1972  at  heights  of  280-370  km  when  the  solar  activity  index  S^q>7 
is  at  levels  exceeding  200. 

To  this  conclusion  may  be  added  three  others,  also  tentative: 

(b)  The  discrepancy,  being  less  towards  the  end  of  the  life  of  1980-43A, 
probably  decreases  as  height  decreases  from  370  to  280  km. 

(c)  The  CIRA  scale  heights  and  densities  are  probably  too  low  over  a 
wider  range  of  heights,  250-500  km,  when  S1Q  7  exceeds  200. 

(d)  The  discrepancy  probably  increases  as  solar  activity  increases. 

These  conclusions  are  also  in  accord  with  the  comparison15  between  CIRA  and  the 
values  of  scale  height  for  early  1958,  obtained  from  Sputnik  2,  Explorer  1  and 
Explorer  3. 

If  the  results  obtained  here  for  1980-81  are  accepted,  the  density  scale 
height  in  CIRA  1972  should  be  about  15Z  higher  for  Sjq  7  *  220  ,  and  this 
amendment  might  be  made  by  taking  a  value  of  exospheric  temperature  substan¬ 
tially  higher  than  recommended  in  CIRA  1972.  A  14Z  increase  in  T—  ,  from  1400  K 
to  1600  K,  Increases  the  density  scale  height  by  9Z,  for  a  height  of  300  km,  so 
an  increase  of  more  than  200  K  seems  to  be  required.  In  CIRA  1972  the  value  of 
T^  is  taken  to  vary  linearly  with  Sjq  7,  but  it  seems  that  the  increase  is 
faster  than  linear  for  Sjq  7  >  200  .  Any  suggested  amendment  can  only  be 
speculative  and  tentative,  but  an  addition  of  100  K  at  T^  ”  1300  K  and  an  addi¬ 
tion  of  200  K  at  T—  ■  1400  K  would  seem  to  be  a  suitably  cautious  revision  that 
might  give  a  closer  approximation  to  reality  during  the  high  solar  activity 
of  1981. 

The  more  recent  MSIS  model16  appears  to  give  higher  densities  than  CIRA 
1972  at  high  solar  activity,  though  comparisons  are  difficult  to  make  because 
the  many  parameters  Involved  are  differently  treated  in  the  two  models  and  the 
tabulations  are  different. 
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7  CONCLUSIONS 

The  orbit  of  NOAA-B  (1980-43A)  has  been  determined  at  40  epochs  between 
October  1980  and  May  1981  from  about  3000  radar  and  optical  observations,  of 
which  the  majority  were  US  Navy  and  N0RAD  observations.  The  average  orbital 
accuracy,  radial  and  cross-track,  was  about  100  m,  with  a  somewhat  better 
accuracy  during  the  final  14  days. 

The  orbital  inclination  exhibits  a  decrease  of  nearly  0.09°  due  to 
atmospheric  rotation,  and  its  variation  was  successfully  analysed  to  derive  two 
good  values  of  atmospheric  rotation  rate,  namely  1.10  ±  0.10  rev/day  for  height 
300  km  (average  local  time  and  season)  and  1.15  ±  0.06  rev/day  for  height  225  km 
(evening  and  average  season).  See  Figs  4  and  6. 

It  has  long  been  known  from  theory  that  the  precession  of  a  satellite's 
orbital  plane  should  be  slightly  altered  by  the  action  of  atmospheric  rotation: 
but  this  is  difficult  to  detect  and  has  not  previously  been  demonstrated.  The 
effect  is  clearly  apparent,  however,  for  1980-43A  during  its  last  8  days  in 
orbit,  and  the  observed  change  agrees  well  with  that  predicted  by  theory  for  an 
atmosphere  rotating  at  1.15  rev/day. 

Ten  values  of  density  scale  height  H  have  been  determined  from  the 
decrease  in  the  satellite's  perigee  height.  These  values,  which  should  be 
accurate  to  about  32,  are  on  average  152  higher  than  indicated  by  the  COSPAR 
International  Reference  Atmosphere  1972  for  heights  of  280-370  km,  although 
values  of  H  obtained  from  previous  orbit  analyses  agreed  well  with  CIRA  1972. 
The  levels  of  solar  activity  prevailing  in  the  years  1980-81  were  higher  than  at 
any  time  since  1958,  and  it  appears  possible  that  CIRA  1972  underestimates  the 
density  scale  height  and  the  density  at  high  levels  of  solar  activity;  very 
little  data  from  1958  is  included  in  CIRA  1972,  so  this  result  is  not 
surprising.  Suggestions  are  made  for  modifying  CIRA  1972. 

The  satellite  NOAA-B  was  regarded  by  NASA  as  a  complete  failure  because  it 
entered  a  wrong  orbit.  Tracking  by  NASA  ceased  soon  after  launch  and  NOAA-B 
never  fulfilled  its  Intended  role  of  supplying  weather  pictures  to  the  world. 

But  it  has  instead  now  given  some  unique  and  useful  information  on  the  rotation 
rate,  density  and  temperature  of  the  upper  atmosphere  at  a  time  when  solar 
activity  was  higher  than  at  any  time  in  the  previous  20  years.  So  NOAA-B  did 
succeed  in  giving  a  picture  of  the  atmosphere,  though  at  higher  levels  than 
intended. 
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ORBITAL  PARAMETERS  FOR  NOAA-B  AT  40  EPOCHS,  WITH  STANDARD  DEVIATIONS 


Table  2  (concluded) 
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Fig  3  Values  of  orbital  decay  rate.  Mo  ,  from  Table  2 
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Fig  5  Calculated  variation  of  inclination  due  to  15th-order  resonance 
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Fig  7  Values  of  n  cleared  of  zonal  harmonic  and  lunisoiar  perturbations, 
with  theoretical  variation  for  A  *  1.15 
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Fin  8  Pengae  height  over  'Earth'  of  radius  6387  km,  after  removal  of  gravitational  perturbations 
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Fifl  9  Valuta  obtained  for  density  scale  height  Hj  at  height  with  CIRA  values  for  comparison 
and  solar  and  geomagnetic  indices  above 
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